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Abstract: Many problems of cold dark matter models such as the cusp problem and the 
missing satellite problem can be alleviated, if galactic halo dark matter particles are ultra- 
light scalar particles and in Bose-Einstein condensate (BEC), thanks to a characteristic 
length scale of the particles. We show that this finite length scale of the dark matter 
can also explain the recently observed common central mass of the Milky Way satellites 
(~ 10 7 Af Q ) independent of their luminosity, if the mass of the dark matter particle is about 
10- 22 eV. 



Keywords: Dark matter, Bose-Einstein condensate, Dwarf galaxies 



© SISSA/ISAS 2009 



http : // jhep . sissa . it/ JOURNAL/ JHEP3 . tar . gz 



Contents 



|]. Introduction |l| 

Minimum mass of galaxies from the minimum length scale |3| 

Minimum length scale of galaxies revisited ^ 

|j. A numerical simulation |6| 

Discussion |8| 

M 

tr^L 

O 

1. Introduction 

Dark matter (DM) is one of the most important puzzles in modern physics and cosmol- 
ogy. Since the presence of DM was inferred from gravitational effects on visible matter 
in galaxies (IJ, a good DM model should explain the structure and evolution of galaxies. 
The dwarf galaxies seem to be the smallest dark matter dominated astronomical objects 
and, hence, are ideal for studying the nature of DM 0], Observational data suggest that a 
typical dwarf galaxy has a size never less than O(10 2 ) pc. This could pose a challenge to the 
structure formation theory based on a simple cold dark matter (CDM) model 0, because 
there is no scale that can be easily accessed observationally in the CDM theory. (However, 
there is surely a scale in the structure formation of baryons in DM halos. Understanding 
theoretically the scale is a major topic of present-day research.) 

Although CDM such as WIMP( Weakly Interacting Massive Particles) with the cosmo- 
logical constant (ACDM) model is popular and remarkably successful in explaining large 
scale structures of the universe, it seems to encounter many other problems such as the 
cusp problem and the missing satellite problem on a scale of galactic or sub-galactic struc- 
tures |J. For example, numerical studies with the ACDM model usually predict a diverging 
central halo DM density pdm(t) °c l/r a ,a ~ 1, where r is a distance from the center of 
a galaxy, while observations indicate that some low surface brightness galaxies have an 
almost flat core density. (However, some dwarf spheroidal satellites of the Milky Way 
have steeper DM density profiles.) The CDM theory also predicts many satellite galaxies 
(O(10 3 )) around the Milky Way, however, only O(10 2 ) satellite galaxies are observed at 
most 

Recently, it was also found that the mass enclosed within a radius of 300 pc in these 
galaxies is approximately constant (~ 1O 7 M0) regardless of their luminosity between 10 3 
and 1O 7 L j8|. This result implies the existence of a minimum mass scale for the dwarf 
galaxies || ||| independent of their baryon matter fraction, in addition to the minimum 
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length scale. This observational fact could be also problematic with the ACDM model, 
which usually predicts many dark matter dominated structures smaller than the dwarf 
galaxies down to 1(T 6 M @. ( Stri gari et al || suggested that this difficulty can be 
overcome if we consider roles of baryonic matter.) 

In this paper, we suggest that a cold dark matter model based on Bose-Einstein con- 
densate (BEC) or scalar field dark matter (SFDM) |l2| can explain this minimum mass 
as well as the minimum length scale observed. Considering observational and theoretical 
uncertainties, we are here concerned with the scales correct within an order of magnitude. 

Let us briefly review the history of the DM theory based on BEC. In 1992, to explain the 
observed galactic rotation curves, Sin [11, ^] suggested that galactic halos are astronomical 
objects in BEC of ultra-light (with mass m ~ O(10 _24 )ey) DM particles such as pseudo 
Nambu-Goldstone boson. In this model the halos are like gigantic atoms, where the ultra- 
light boson DM particles are condensated in a single macroscopic quantum state ip. The 
quantum mechanical uncertainty principle prevents halos from a self-gravitational collapse. 

Lee and Koh [14, [T^] suggested that this condensated halo DM can be described as a 
coherent scalar field JTI| (dubbed as the boson halo model or, more generally, SFDM model 
later [16| ). On the other hand, in usual CDM models, the wave functions of DM particles 
do not overlap much and the particles move incoherently. Thus, the key difference between 
our model and the usual CDM models is the state of DM particles rather than DM itself. 
Later, similar ideas were suggested by many authors in terms of the fuzzy DM, the fluid 



18L H9 20, 21, 22, 23, 24, 25, Efl, 27 



DM, or the repulsive DM [ 
(See p] for a review.) Many related models like monopole failing models 
massless scalar field models [FT; 



, 31 



351, unstable 



and quintessence models |19|] are suggested too. There are 
also extensive literature on, for example, the collision of BEC-dark matter structures [36, 
|37|1 , non-spherical collapse |3g] and virialization processes [pSfl , and the collapse of SFDM 
fluctuations in an expanding universe [^, 39, 40, 41, |4^[. (For a review sec |i"5| , |i~7| , 28] and 
references therein.) 

In these models (BEC/SFDM model hereafter) the formation of DM structures smaller 
than the Compton wavelength (0(1) pc) of DM particles is suppressed by the quantum 
uncertainty principle. (We will see that the actual length scale is somewhat larger.) It was 
shown that this property could alleviate the problems of the ACDM model [23, |25|, 



such as the cusp problem [45] and missing satellite problem. It is also suggested that this 
model can well explain the observed rotation curves 
and the size evolution of the massive galaxies [f 



], collisions of galaxy clusters ||37[| 
Note that the BEC/SFDM particles 



behave like CDM [48] at the scale larger than galaxies during the cosmological structure 
formation due to a small velocity dispersion of the Bose-Einstein condensated state. Thus, 
this model safely satisfies the criteria of the well established structure formation theory 
above the galactic scales. 

In this paper, we show that the minimum mass of galaxies can be explained if DM 
is in BEC. In Sec. 2 we show the relation between the minimum mass and the minimum 
length scale of galaxies. In Sec. 3 we review the proposed origin of the minimum length 
scale in the BEC/SFDM model. In Sec. 4 we present a result of a simple numerical study 
supporting our theory. Section 5 contains discussion. 
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2. Minimum mass of galaxies from the minimum length scale 



The characteristic length scale (0(1) ~ O(10 3 ) pc) for the clustering of dark matter is the 
key feature of the BEC/SFDM |0, |l2| and its variants. It is very interesting that this 
minimum length scale of DM halo seems to be already observed |^j . There are no known 
stable galaxies with a half-light radius smaller than 120 pc while the maximum size for star 
clusters deprived of DM is about 30 pc. Furthermore, observational data indicate there is 
a distinction in phase-space density between a star-cluster and a dwarf galaxy [^9|. Since 
galaxies, especially dwarf galaxies, are highly DM dominated objects, it is plausible that 
the minimum mass of galaxies observed is connected with DM rather than visible matter 
in the galaxies. Thus, it is natural to think that the minimum mass of galaxies is somehow 
related to this finite length scale. 

i r~ i 

First, let us briefly review the BEC/SFDM model, in which a galactic dark halo is 
described by a quantum wave function vp(r) (or scalar field) of the non-linear Schrodinger 
equation with the Newtonian gravity 0]; 

ft 2 

ihd t ip = Eip = -_vV + V(rWr), (2.1) 
2m 

where E is energy of a DM particle and r is the distance from the halo center. The 
gravitational potential of the halo is given by 

f r Cm r r> 

V(r)= dr'— dr"4nr" 2 (p vis (r")+p DM (r")) + V , (2.2) 
Jo r Jo 

where DM density pdm(^) = Mo\ip(r)\ 2 and p v i s is the visible matter (i. e., stars and gas) 
density. Mq is a mass parameter and Vq is a constant to make ^(oo) = 0. These equations 
can be also derived in the SFDM theory in the Newtonian limit. For simplicity, we assume 
there is no self-interaction between the scalar particles except for the gravity. 

Let us find the size of a stable configuration of the halo. From Eq. ( |2.1| ) the energy E 
is approximately given by 

B (0-o^2+/ dr'^ dr"^r" 2 (p ms (r")+PDM(r")l (2.3) 
zml; J r j 

as a function of a halo length scale £. A stable ground state configuration of the halo suitable 
for dwarf galaxies (a zero-node solution) can be approximately found by extremizing the 
energy by £ @; 

, ^ i h 2 GMm , N 

dE{g) dt ~ -3 + — — = 0. 2.4 

m^ 6 t^ z 

Here, 

M = f dr"^r" 2 {p ms {r")+p DM {r")) (2.5) 



o 



is the total mass within £ of the galaxy consisting both of the DM and the visible matter. 



The condition in Eq. (2.4) satisfies at |ll 



h 2 c 2 \ 2 c 

^ = GMm 2 = 4vr 2 GM' ^ 
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where c is the light velocity and A c = 2irh/mc is the Compton wavelength of the particles. 
From this equation one can obtain the mass within £, 

M (0 * T^r~2- ( 2 - 7 ) 

G£m z 

It is well known that, in the BEC/SFDM theory, the size of the stable configuration 
£ could not be arbitrary small. Let us denote this minimum value of £ as £ c which is at 
least A c . Then, there appears a minimum mass for the DM halo, M c = M(£ c ) = h 2 /G^m 2 
corresponding to £ c . If £ c = A c ~ 1/m and m ~ 10~ 24 eV, the minimum mass for the DM 
halo becomes M c = hc/2TrGm ~ 1O 13 M , which is too large for dwarf galaxies. However, 
note that the free parameter mass m ~ 10~ 24 eV is just one of rough estimates. In the 
literature on the BEC/SFDM, m was usually suggested to be in the range 0(W~ 2e )eV ~ 
0(W~ 22 )eV to solve the problems of CDM models. 

Conversely, from the observed values of M c and £ c one can obtain approximate mass 
of DM particles 



h 2 



m 



From the observed value M c ~ 10 7 M Q , one can obtain m ~ 5.4 x 10~ 22 eV for £ c = 
300 pc from Eq. (|2. §|), For this m, A c ~ 0.075pc. (In the next section we will discuss 
why £ c is larger than A c .) Very interestingly, this m value is similar to the value m = 
10~ 22 eV ~ W~ 23 eV required to solve the cusp problem and to suppress the small-scale 



power |22|, 26 , 511] . Thus, m ~ 10 22 eV can solve all known small scale problems mentioned 



above, which was discussed by Hu, Barkana and Gruzinov in Ref. p^j . 
3. Minimum length scale of galaxies revisited 



In fact, even in the original work [11, 12] it had been already noted that the typical length 
scale (0(kpc)) of galaxies is somewhat larger than the Compton length A c . This is related 
to the fact that the halos are basically non-relativistic objects. £ c could be comparable 
to A c only when astronomical objects in consideration are extremely relativistic like black 



holes. It is clearly not the case for galactic halos. £ in our equation (2.6) was identified 
as the gravitational Bohr radius ~ kpc in [|llj. In |[^] it was argued that de Broglie 
wavelength A^b, rather than the Compton length, is more adequate for the typical size of 
DM structures. (Similarly, in condensed matter BEC systems the coherence length or the 
healing length, usually comparable to A^b, rather than A c determines the spatial size of 
BEC density fluctuations |52[] .) 

There are several ways to obtain £ c slightly larger than A c . The first and most plausible 
one is the quantum Jeans scale rj ~ 0({Gpm 2 )~ l l 4: ) |5^, 22, ^0| for SFDM, which is the 



geometric mean between the virial dynamical scale and the Compton scale. This scale can 
be derived from the cosmological evolution equation for the SFDM density perturbation 5 
with a wave vector k [22, 54 1 

d 2 5 + 2Hd t S - (akG P - J^-j) 5 = 0, (3.1) 
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where H is the Hubble parameter, p is the mean DM density and a is the scale factor. 
The second term in the parenthesis represents a contribution from quantum pressure of 
the DM. The density perturbations are stable below rj and behave as ordinary CDM 
above this scale. This fact makes the BEC/SFDM an ideal alternative to the CDM. 
rj determines the minimum length scale during the cosmic structure formation, which 
could be larger than A c . (This quantum Jeans instability is also known in the studies of 
boson stars [55, According to the standard structure formation theory, a physically 

meaningful value of rj can be fixed at the matter-radiation equality having a scale factor 
a eq ~ 1/3200. At this time p ~ 3 x 10 10 M Q a~q /Mpc 3 []57|], and the quantum Jeans scale rj 
above becomes O(10 2 )pc for m — W~ 22 eV. Hence, £ c ~ rj is plausible. Second, £ c can be 
the coherence length sj Mp /ipah/mc determined by the DM scalar field value ipo at the halo 
center [^], which might be related to some symmetry breaking. For a simple quadratic 
potential p ~ m 2 ^ 2 ,, the coherence length again becomes of order rj |E(J. This implies 
V>o ~ 10 _8 Mp, where Mp is the Planck mass. Third, if there is a self-interaction term Xcj) 4 
between DM particles, a new length scale £ c = 0(A 1 / 2 A c M P /m) emerges |17|]. Finally, one 
may also consider a thermal de Broglie wavelength for £ c [j5Q(] - Thus, £ c somewhat larger 
than A c is not new in these models and theoretically possible. 

We assume one of these scales plays a role of £ c , though the quantum Jeans scale seems 
to be most plausible for our purpose. Note that all these scales are determined not by the 
properties of individual halos but by those of the DM particles such as mass, coupling and 
mean density of the DM particles. Therefore, ^ c ~ O(10 2 ) pc > A c could be a universal 
quantity H. In short, the nature of the DM particles fixes £ c , which in turn decides M c . 

To make a DM dominated structure, the total mass of the structure should be large 
enough to gravitationally attract DM particles having intrinsic momentum of 0(h/^ c ) due 
to the uncertainty principle. The critical mass M c is just the minimum value to do this. 
From the virial condition, i. e., by equating kinetic energy ((h/^ c ) 2 /2m) with potential 
energy (V = GMm/^ c ), one can check that M c in Eq. ( |2,7| ) is of correct order. It is 
important to note that, since £ c is the property of DM itself, M c should be independent of 
the fraction of visible matter, hence the luminosity, as long as the fraction is small. This 
seems to be the basic physics behind what was observed. 

The thermal history of this dark matter is also considered by several authors [58, 50, 59]. 
Since the dark matter particle is ultra-light, its thermal velocity could be highly relativistic, 
and some concerns about the BEC transition may arise. There are basically two options. 
First, BEC/SFDM might simply have been generated by a non-thermal process like axions. 
Second, even when the DM particles are in thermal equilibrium, contrary to common 
wisdom, a BEC could form for relativistic bosons [SO, [H], |62|, 59]. In |5£| it is shown that 
a cosmological BEC always exists for the relativistic SFDM, if m < 10~ 14 eV |30], p2j| . 

After the BEC phase transition at the temperature T c , almost all DM particles are 
in a ground state and move coherently rather than randomly. The BEC ground state 
are favored against thermally excited states owing to the Bose-Einstein statistics of boson 
particles p0[ . According to [ |59| ] the present DM number density n ~ 10 15 ey 3 is much 
larger than the maximum charge density allowed to excited states; mTg/S ~ 10 _30 ey 3 , 
even when the present temperature of the DM particles is as high as Tq ~ 10~ 4 eF for 
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m ~ 5.4 x 10~ 22 eV. (See Eq. (13) of the reference.) Hence, almost all BEC/SFDM 
particles are in the ground state now even in this case. 

Once the relativistic BEC is formed, SFDM particles in the ground state have only a 
small quantum velocity dispersion ~ h/m£ and could be bounded in a self-gravitational 
potential (if M > M c ) even though the temperature of the condensate could very high 
compared to the mass of the particles. The quantum velocity dispersion Av ~ h/m^ c ~ 
'WCc ~ 10~ 4 c ~ 0(W)km/s is smaller than the escape velocities of dwarf galaxies. For 
stability analysis, we can think of the boson halos as boson stars [12, |l7]]. Their stability 



with respect to various kinds of perturbations or gravitational collapse is shown in many 
works Q. All these facts indicate the stability of the dark halos with the BEC/SFDM Q. 
There are also universal profiles for boson stars, which are attractors of the collapse of quite 

EE 



arbitrary initial density fluctuations [64]. 

4. A numerical simulation 

To be more concrete, we perform a numerical study using the shooting method |l2j for a 
toy DM halo. Similar numerical work for cases without visible matter has been extensively 



performed in the studies of boson stars [16|. The dimensionless form of Eq. ( |2.1[ ) and Eq. 
( p.2[ ) can be written as 

f V?V = {a 2 + p ms ) 

[V 2 P a = 2(V - E)a ' 1 1 

o 



Here we have introduced dimensionless quantities |12j 

a = v / 47Ge-^V/c 2 , (4.2) 
f = mcr/h, 
i = mc 2 t/h 

so that all other barred quantities are dimensionless. For example, in this unit, £ = mc^/h, 
V = V/c 2 and E = E/mc 2 . 

For visible matter of dwarf galaxies, we use the dimensionless version of the empirical 
density in Ref. pf [ 

/ - \ 3/3 

p vis = 3/3M vis ( ^— , (4-3) 
\r c + rj Qnr^lr + r c ) 

where (3 = 2 and M v i s is a parameter proportional to the visible matter fraction. It has been 
shown that this p v i s with the BEC /SFDM can successfully reproduce the observed rotation 
curves of dwarf galaxies |46[| . We choose the visible core size f c = 664 corresponding to the 
physical size r c ~ 40 pc. Considering the observational data, we choose £ c = 3985 (£ c = 
300 pc). One can calculate a dimensionless total mass within £ c , 

M tot = / A7rr 2 (a(r) 2 +p vis (f))dr, (4.4) 
Jo 

which corresponds to a physical mass 
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Figure 1: (Color online) The dark matter density a 2 (blue dotted lines), the visible matter density 
Pvis (green thick lines) and the gravitational potential V (red dashed lines, in units of c 2 , scaled 
by 3 x 1CP 8 ) for a model dwarf galaxy as functions of distance f from the halo center. Despite 
of wide variations of a visible matter fraction, the 3 profiles have a universal core size r c = 4000 
(corresponding to the physical size r c ~ 300pc) and the total mass within this size Mtot — 0.00019 
(corresponding to 4.75 x 1O 7 M ) is similar for all 3 cases and the gravitational potential functions 
almost overlap. Here, f — 1 corresponds to a physical distance 1/m ~ 0.075 pc and a = 1 to the 
field value \tp\ = l/V^brG. 

Fig. 1 shows the result of our numerical study with boundary conditions dV /df(0) = 
0, V(oo) = and da/df(0) = 0. For other parameters, we consider 3 cases with the 
parameters (M vis = 5 x 10~ 6 ,o-(0) = 3.6 x 10~ 8 ,V(0) = -2.18 x 10" 8 ,^ = -1.7 x 
l(r 9 ),(M» s = 3 x 10- 5 ,a(0) = 3.4 x 10~ 8 ,V(0) = -2.12 x lO" 8 ,^ = -1.3 x 1(T 9 ), and 
(M vis = 1 x lfr 4 ,cr(0) = 2 x irr 8 ,F(0) = -2.097 x 10~ 8 ,^ = -1.28 x l(r 9 ), respectively, 
from the top to the bottom for <r(0) (reversely for p v i s ). ( One can easily recover the 
dimensionful quantities by using the transformation in Eq. (|4.2| ) and below. For example, 
E = —1.7 x 10 -9 corresponds to the real energy E = Emc 2 = —9.18 x 10 -31 eV.) Three 
profiles of the potential functions almost overlap as our theory expected. For a similar 
mass M t ot — 0.00019, the shapes of three V(r) are similar, regardless of M v i s , as long as 
DM is dominant. The converse is also true. 

The result confirms our arguments below Eq. Q2.7| ). Despite of wide variations of visible 
matter fraction, the profiles of DM and potential well have a universal core size (£ c ~ 300 pc) 
and the total masses within 300pc, M tot = (0.000189, 0.000190, 0.000189), are similar for all 
3 cases. (In the physical units these are M tot ~ 0.00019MJ, /m ~ 4.75 x 10 7 M o ~ M c ). The 
masses of the visible matter are (1.25 x 10 6 , 7.5 x 10 6 , 2.5 x 1O 7 )M , respectively. Thus, the 
densities of dark and visible matter are similar to the observed values. In our simulation, 
what really matter were M c and £ c not the composition of matter within £ c . This confirms 
that there is a universal minimum mass of galaxies composed of the DM and ordinary 
matter corresponding to the characteristic length scale in the BEC/SFDM theory. This 
consistency of the theory and observations is non-trivial because the numerical solutions 
of Eq. ( fl.ip are very sensitive to the boundary conditions. On the contrary, the results are 
not sensitive to the form of p v i s . We have checked that a similar conclusion can be derived 
for p v i s with (3=1. When the DM is subdominant, the universal minimum mass and size 
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are not shown clearly. This is also expected because the universality is a property of the 
DM not that of visible matter. 



In [64] it is shown that under isolation conditions the rotation curve would be Keplerian, 
unless an appropriate spatial scale is chosen. We have used the cut-off of the scalar field 
profile at the boundary f = f o = 22581 (corresponding to about 1700 pc) so that the density 
of the field equals the present average DM density of the universe po ~ 3 x 10 10 Mq/Mpc 3 , 
which corresponds to cr(ro) — 10~ n in the dimensionless unit. This is necessary, because 
the DM density of galaxies is usually non-zero even for large r. With ip ~ 10 _8 Mp and 
m ~ \0~ 22 eV, we can successfully reproduce the observed halo DM density m 2 if) 2 ~ 
10" 4 ey 4 ~ 10 7 M o /(100pc) 3 for dwarf galaxies as well as £ c ~ O(10 2 pc). Since above the 



galactic scale, SFDM basically behaves as CDM 22], we think that the usual bottom-up 
hierarchial merging process happened for clusters of galaxies. 

5. Discussion 

rTj 

We have shown analytically and numerically that the BEC/SFDM theory can explain the 
three observed properties of dwarf galaxies, i.e., the minimum length scale, the minimum 
mass scale, and their independence from the brightness. On the other hand, it is difficult 
to explain all these properties in a single scenario in the CDM context, although there are 
proposed solutions to some of these problems f66| , |6^| relying on roles of baryon matter. 
Another merit of our approach is that the stable DM configurations are not much dependent 
on the complicated galaxy merging history, because the equilibrium condition forces the 
halos to rearrange their DM distribution so that they have a universal form regardless of 
their history. This could explain the observed universality of the mass and density profile 
independent of possible merging history. 

In conclusion, the BEC/SFDM could be a compelling alternative to the standard 
CDM, because it could not only solve some known problems of CDM model such as the 
cusp problem and the satellite problem, but also has a possibility to explain recent obser- 
vational mysteries of galaxy evolution in a simple way, thanks to its wave nature and the 
characteristic length scale. Further tests and future observations may reveal whether CDM 
is in a BEC. 



This work was supported by the topical research program (2009-T-l) of Asia Pacific Center 
for Theoretical Physics. 
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